Granular Dynamics during Impact 
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We study the impact of a projectile onto a bed of 3 mm grains immersed in an index-matched 
fluid. Specifically, we vary the amount of prestrain on the sample, strengthening the force chains 
within the system. We find this affects only the prefactor of linear depth-dependent term in the 
stopping force. Using a laser sheet scanning technique to visualize internal grain motion, a high- 
speed camera, and particle tracking, we can measure the trajectory of each grain throughout an 
impact event. Microscopically, our results indicate that weaker initial force chains result in more 
irreversible, plastic rearrangements during impact, suggesting friction does play a role. However, 
by comparing results of wet and dry systems, we infer that friction is not the primary cause of the 
linear depth-dependent term. These results suggest the depth-dependent mode of energy dissipation 
is likely acoustic in nature, governed by enhanced normal contacts between grains. 



Take a run on the beach, and your foot strikes a gran- 
ular material in much the same way an asteroid hits a 
planet. Even though this interaction is commonplace, 
the physics of it remain largely mysterious. Previous 
work has mostly focused on the intruder's dynamics and 
has focused on continuum models, offering little insight 
into the microscopies, where the actual grain rearrang- 
ments take place [1-9]. It's not surprising the microscop- 
ies are ignored; granular materials are difficult to study. 
In addition to being opaque, granular materials are not 
well-behaved: they are heterogeneous, can behave either 
like liquids or solids, and exhibit shear localization and 
jamming phenomena [10]. To make predictions on prac- 
tical applications of granular impact, such as how far an 
asteroid will penetrate into soil, it is clear we first need 
a more complete understanding of the entire system. 

When an intruder goes into a granular material, the 
material behaves somewhat counterintuitively: it be- 
comes stronger when it is struck harder. The stopping 
force increases as the impact energy (in this case, drop 
height) increases, resulting in a shorter stopping time for 
higher energy impacts. One consequence of this behavior 
is a a typical scaling law found: 



d ~ H l ' z 



(1) 



, where d is the penetration depth and H is the total 
drop height, the initial height above the bed plus the pen- 
etration depth [8] . As the impact energy (drop height) is 
increased, it is dissipated over a (relatively) shorter pen- 
tetration depth. However, this scaling does not appear 
to be universal - different studies report different scaling 
exponents for similar systems, underscoring a need to see 
into the details [8]. Further, much work has focused on 
the same initial state: loose granular matter, stirred or 
fluffed by air flow. More recent work has shown that con- 
trolled modifications to the packing fraction, cohesion, 
and interstitial gas can produce drastically different out- 
comes [11-13]. 



In the loosely packed system, Katsuragi and Durian 
[1] interpreted the depth scaling as the result of an 
empirically-determined force law: 



EF 



-mg + k\z\ + mv 2 jd\ 



(2) 



, where k and d\ are constants obtained by fits. The 
mv 2 jd\ term, is an inertial drag term, arises from mo- 
mentum transfer during grain-intruder and grain-grain 
collisions. The form of this force law is generally agreed 
upon for shallow impacts, though the friction-like k\z\ 
term saturates with deeper impacts [7]. Katsuragi and 
Durian describe the k\z\ term as Coulomb friction, but 
its magnitude is 20 times higher than one would expect 
from Coulomb friction. This raises the question whether 
Coulomb friction is even responsible for this term. On 
the other hand, Hou and coworkers [5] found the depth- 
dependent term is the correct order of magnitude consis- 
tent with grain-grain friction. Similarly, Tsimiring and 
Volfson [3] found that friction was responsible for energy 
dissipation in their simulated system. However, in an- 
other simulation, Seguin et al [14] found that the dynam- 
ics of impact can be reproduced with no friction between 
grains. 

To settle some of these debates, we must look into 
how the material dissipates the impact energy at grain- 
grain contacts by analyzing rearrangements at the scale 
of the grain motion. Our hypothesis is that grain level 
"failure" or rearrangements will be strongly dependent 
on the initial state, or "fabric" of the granular material. 
Granular materials have mesoscopic features that signify 
their fabric: force chains. These are filamentary networks 
of contacting grains that bear the weight of the system, 
transmit energy along their length, and resist shear and 
buckling due to the friction between grains [15]. Thus by 
systematically varying the initial force chains within the 
system, we can better understand the role of this network 
in the dynamics of impact. 

In this Letter, we systematically vary the initial state 
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FIG. 1. A schematic of the experimental setup. The beads 
are in an index- matched fluid (DMSO) with a fluorescent dye, 
allowing illumination of the sample in the impact plane. The 
sample is strained via displacement of the back wall. 



by pre-straining the sample, in order to vary the strength 
of force chains within the sample. We analyze how this af- 
fects the macroscopic scaling of the intruder's dynamics, 
and thus the stopping force. We investigate the micro- 
scopic behavior of the granular material. To do this, we 
look inside a granular bed during impact, tracking par- 
ticles near the intruder. This allows us to characterize 
mesoscopic failure events during the impact. 

Our experimental setup is shown in Fig. 1. Our grains 
are 3 mm glass borosilicate spheres (Glen Mills). They 
are poured to fill a clear box, approximately 15 cm on 
each side. The box has one moveable sidewall to apply 
strain to the system. The box is then filled with a mixture 
of DMSO, water, and Nile Blue 690 perchlorate dye. The 
mixture is tuned to match the index of refraction of the 
grains, so that the index mismatch is less than 0.005. 
A laser sheet is set to illuminate the plane of impact, 
resulting in bright fluid and dark grains [16]. 

We image this slice during impact using a high speed 
PCO.edge camera. The laser sheet and moveable sidewall 
are positioned by stepper motors. The impactor is a 1 
inch carbon steel sphere, released via an electromagnet. 

With the moveable sidewall, we can build up force 
chains in the sample by (very slightly) compressing the 
material from the side. The force chains will then be 
biased perpendicular to the direction of impact. As we 
start at a loosely packed state (<j) « 0.55), we expect slight 
compaction when straining the sample. We also keep the 
height of the sample fixed with a top plate. Thus max- 
imum volume fraction increase is approximately 0.003, 
less than the magnitudes in [13] needed to have an ef- 
fect on dynamics. So with higher values of the strain, we 
create a stronger force chain network only, not a denser 
packing. After the sample is strained, the projectile is 
dropped into the bed and the impact is recorded. We 
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FIG. 2. Depth vs total height for many different sample 
preparations. The fits shown are to power laws, with ad- 
justable prefactor and exponent. Each point corresponds to 
at least three experimental runs. The dry and wet unstrained 
samples show the same exponent ~0.4, and differ in prefactor 
by about 15%. Prestraining the sample causes some deviation 
from the power law at low drop height. A slight increase in 
the scaling exponent is also seen, up to «0.5. 



restore grains to an unstrained state by reversing the 
compression, and stirring the sample for several minutes. 

From the videos we can extract 1) the initial and final 
position of the intruder, 2) the instantaneous position 
(and velocity) of the intruder, and 3) the positions (and 
velocities) of the grains, using our established particle 
tracking routines [17]. 

We examine the scaling of d, the penetration depth, 
with the total drop height H, in Fig. 2. Each data point 
corresponds to five experiments. For our unstrained sys- 
tem (dark blue circles), we see that indeed a power law 
seems to match the data quite well. The best exponent 
for this scaling is ~0.4, a departure from the simple 1/3 
scaling but well within the range of exponents seen for 
various shallow impacts (^0.25-0.50). 

We have also looked at our system and compared it 
to a dry system (light blue squares). Surprisingly, the 
depth scaling does not change for wet vs dry, the only 
difference is in the prefactor, about a 15 percent differ- 
ence. The fluid dynamics are the same, our fluid is low 
viscosity and the impact speeds are fast, meaning that 
flow in the channels between particles is well out of the 
viscous regime. In fact the Reynolds numbers for the 
dry (air) and wet (DMSO) system are quite similar, and 
both greater than 1. Second, We are also in a fully satu- 
rated state, so surface tension effects are negligible [12]. 
In fact, the decreased impact into the wet material can 
be explained solely in terms of the increased mass (about 
30 percent increase) due to the presence of the fluid, as 
will be detailed later. 

We find that straining the sample greatly changes the 
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FIG. 3. Analysis of the intruder trajectory, (a) Three sample 
trajectories for a fixed drop height with differing amounts of 
prestrain. From the trajectories, we are able to find a value 
of (b) d\ for each strain value and (c) k for each strain value. 
The dashed and dotted lines show the values of d\ and k for 
the unstrained wet impact, if the additional mass of the fluid 
is scaled out. 



power law scaling of depth vs. drop height (Fig. 2). Even 
small strain results in a much shallower impact. Further, 
the exponent is increased, though only to a maximum 
value of ^0.5. This suggests that with stronger force 
chains, the impact hardening of the material is less, per- 
haps because the material is already somewhat hardened. 
We measure the trajectories of the intruders, an exam- 
ple is given in Fig. 3a for three different strains. From 
this we can analyze the stopping force. We rewrite equa- 
tion 1 in a more general way: 



SF 



-mg + h{z) + f 2 (z)z 2 



(3) 



By recasting this equation in terms of kinetic energy, 
we are able to obtain exact results for /i and f 2l using 
only the position data and its first derivative. The details 
of this procedure may be found in [18]. For our data, we 
find that Eq. 1 is suitable to describe the motion at 
all strains: fi(z) ~ z and f2(z) = constant. The slight 



change in depth scaling exponent with strain should in 
principle slightly alter the form of one or both of these 
terms, but our resolution is not sufficient to measure this 
small effect. 

Thus, we calculate the prefactors k and d\ and see how 
they depend on strain, as shown in Figs. 3b and 3c. Each 
point is calculated using data from several drop heights. 
We see the inertial constant, di, does not vary with in- 
creasing strain. On the other hand, the "frictional" con- 
stant k does monotonically increase with strain. Our 
value of d\ is factor of 2 smaller than measured in [1], 
perhaps due to the difference in bead/intruder size ratio, 
as well as the presence of the fluid. Our (unstrained) 
value of k is in agreement with the value found in [1]. 

Also included in Fig. 3b and 3c is data for the dry, 
unstrained sample, k should necessarily scale with the 
mass surrounding the intruder [1, 5] due to its hydrostatic 
pressure-like form and d\ should scale inversely; more 
mass means more material must be moved [1]. When 
properly rescaled with the added weight of the wet sys- 
tem, we find that k and d\ are independent of fluid im- 
mersion, as shown in Fig. 3b and 3c. Further, by rescal- 
ing the wet value for k and d\ in this fashion, and numer- 
ically integrating Eq. 2, we reproduce the depth scaling 
for the dry result. Thus, the difference in /c, di, and scal- 
ing prefactor, are primarily the result of the increased 
mass of the system. 

We have shown that prestrain, which reorients and 
likely enhances the force chains between particles leads 
to an increased kz term in the stopping force. When en- 
hancing the force chains between particles, we enhance 
two forces. First, the normal force at particle contacts is 
enhanced. And because of this enhanced normal force, 
the tangential friction force between particles is also in- 
creased, creating force chains that are harder to buckle. 
It is unclear whether the increase in k is 1) due to a 
strengthened elastic response of the material (higher nor- 
mal forces create an effectively higher modulus [19] or 2) 
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FIG. 4. Particle tracks of the grains during impact. The 
tracks go from red (early) to blue (late). Grains below the 
intruder move down, and grains to the side of the intruder 
move out and up. The frame-to-frame motion is not always 
smooth, due to the frustrated nature of the packing. 
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to the side of the impactor. This flow field has been 
characterized before only in 2D systems, but we find it 
in good qualitative agreement with them [20, 21]. 

We then measure the signatures of buckling force 
chains [15]: local plastic rearrangements of particles rela- 
tive to their neighbors, using the quantity £^ in [22-24]: 



D 2 



liYlAdiM-Eidij]}* 



(4) 



FIG. 5. Characterizing the nonafflne motion, (a) Spatial map 
of D^iri for three different strains at a fixed drop height. The 
amount of nonaffine motion near the intruder decreases as the 
strain is increased, (b) For three drop heights and five strains, 
the average D^iri value is normalized by the kinetic energy 
of the intruder just before impact. The relative amount of 
nonaffine motion is the greatest for lowest drop heights, where 
the kz term is most important. 



due to a weaker plastic response of the material (higher 
friction creates force chains that less easily buckle). In 
short (ignoring the inertial term, which is less important 
for shallow impact), energy can dissipate two ways, by 
acoustically traveling through the normal contacts, or by 
frictional dissipation at the contacts. Both modes can co- 
exist; and we cannot measure the acoustic signal in our 
system. There is new evidence [20] that it plays a role. 

However, the friction between beads might also be in- 
fluenced by the presence of the interstitial fluid. We 
have performed two measurements of friction with our 
beads. We find similar static coefficients of friction 
i^dry = 0.36 ± 0.02 and fi wet = 0.39 ± 0.02 by testing 
the angle of repose. But a test of the sliding friction, by 
dragging two beaded plates across each other, tells a dif- 
ferent story: jidry — 0-35 ± 0.02 and ji wet = 0.26 ± 0.02. 
This dynamic friction is more important for the impact 
environment, but seems to have little effect on the dy- 
namics. Further, this is more of a measure of the mi- 
croscopic particle-particle friction in the system. If the 
friction influences the kz term at all, then we would ex- 
pect k (adjusted for the mass of the fluid) to be smaller 
for the wet case, and we find them to be the same (see 
Fig. 3c, dotted line). This present a strong case that 
k does not primarily represent particle friction. Instead, 
we propose elastic dissipation is primarily responsible for 
the kz term, since the modes travel through the grains 
themselves and are insensitive to the fluid. Further, sim- 
ilarly to [1], our value for k is about 20 times higher than 
one would expect if friction were the sole cause. 

By looking at the trajectories of the particles, we can 
gain further insight into energy dissipation in the system. 
We track the motion of the particles in the plane. An 
example tracked data set is shown in Fig. 4. We see 
the trajectories are characterized by downward motion 
below the impactor, and upwards and outwards motion 



-^min i quantifies the nonaffine motion of j particles in 
the neighborhood around a given particle i after remov- 
ing the averaged linear response to the strain, given by 
tensor Ei\ a smaller £>^ in indicates more affine motion. 
The vector dij is the relative position of i and j, Adij is 
the relative displacement after a delay time At, set to be 
1 frame in this case. 

We find a high E^ in i zone near the projectile (Fig. 
5a); this result is in agreement with 2D simulations of 
impact in [21]. By looking at the average value of D^^ i 
in our field of view, we can compare how the plastic rear- 
rangements (nonaffine motion) varies with initial sample 
strain (Fig. 5b). We find that increasing the initial strain 
results in decreased nonaffine motion. This suggests that 
the buckling of force chains decreases as the chains get 
stronger, suggesting friction is a factor in this system. 
Indeed, in [21], Kondic et al found that by introducing 
friction in their simulation, the nonaffine motion was de- 
creased. 

In Fig. 5b we also show that £>^ in is most disturbed by 
strain for low-energy impacts. This lends further support 
to the notion that friction is important - at lower impact 
impact energies the kz term becomes more important. 
By increasing /c, the plastic response signifying network 
breakage decreases. 

In summary, we have shown that the initial prepara- 
tion of a granular system strongly affects its failure under 
impact. This is the first study of uniaxial strain on such 
as system. We have characterized how strain affects the 
stopping force in the material, showing an increase in 
the so-called frictional term. By measuring bead-scale 
motions, we see that friction does play a role; the trends 
in nonaffine motion support this. 

However, our results shows that this depth-dependent 
term is primarily to particle-particle friction. While 
particle-particle friction is difffernt for wet and dry sys- 
tems, this term does not change. Instead, acoustic dissi- 
pation due to enhanced normal contacts appears to be a 
factor [20]; and these normal contacts will increase with 
depth and strain in the same way the frictional contacts 
do [15]. Perhaps then, for a system of hard spheres, 
normal modes of dissipation are primary and tangential 
modes are secondary, conspiring to create the kz term. 
This aligns with intuition: these spheres are high modu- 
lus and low friction. The exact form of f\ (z) may depend 
on the relative importance of these forces in the material; 



tuning the frictional or acoustic properties of the grains 
should yield further insight into the exact nature of this 
term. 
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